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Abstract 

The production of D"^ and D'^ mesons has been measured with the ZEUS detec- 
tor at HERA using an integrated luminosity of 133.6 pb~^. The measurements 
cover the kinematic range 5 < < 1000 GeV^, 0.02 < y < 0.7, 1.5 < < 15 
GeV and \ri^\ < 1.6. Combinatorial background to the D meson signals is 
reduced by using the ZEUS microvertex detector to reconstruct displaced sec- 
ondary vertices. Production cross sections are compared with the predictions of 
next-to-leading-order QCD which is found to describe the data well. Measure- 
ments are extrapolated to the full kinematic phase space in order to obtain the 
open-charm contribution, Fg'^, to the proton structure function, F2. 
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1 Introduction 



Charm quarks are copiously produced in deep inelastic scattering (DIS) at HERA. At 
sufficiently high photon virtuality, Q^, the production of charm quarks constitutes up to 
30% of the ep cross section [1,2]. Previous measurements of D*^ cross sections [1-5] 
indicate that the production of charm quarks in DIS in the range 1 < < 1000 GeV^ is 
consistent with the calculations of perturbative Quantum Chromodynamics (pQCD) in 
which charm is predominantly produced via boson-gluon fusion (BGF). This implies that 
the charm cross section is directly sensitive to the gluon density in the proton. 

A charm quark in the final state can be identified by the presence of a corresponding 
charmed hadron. In this paper a study of the production of two such charmed particles, 
the and D^/D^ mesons, is presented. The mesons are reconstructed using the decays 
D'^ K~7T^7T^ and K~'it~^, which are chosen as both contain charged particle^ 

which are well reconstructed in the ZEUS detector. The proper decay lengths are of the 
order 300 /im and 100 /im for the and D^, respectively and can be measured [6, 7] 
with appropriate silicon trackers such as those at HI and ZEUS. 

Measurements of the and cross sections are presented with improved precision and 
in a kinematic region extending to lower transverse momentum, p^, than the previous 
ZEUS results [8] ; this is made possible through the use of the precision tracking provided 
by the ZEUS microvertex detector (MVD). Single-differential cross sections have been 
measured as a function of Q^, the Bjorken scaling variable, x, p^, and the pseudorapidity, 
1]^, of the D mesons. The cross sections are compared to the predictions of a next-to- 
leading-order (NLO) QCD calculation using parameterisations of the parton densities in 
the proton which were determined from fits to inclusive DIS measurements from ZEUS 
and fixed-target experiments. The cross-section measurements are used to extract the 
open-charm contribution, Fg^'^, to the proton structure function, F2. 



2 Experimental set-up 

The analysis was performed with data taken from 2004 to 2005 when HERA collided 
electrons with energy = 27.5 GeV with protons of energy Ep = 920 GeV. The results 
are based on an e~p sample corresponding to an integrated luminosity of 133.6 ± 3.5 pb~^. 

A detailed description of the ZEUS detector can be found elsewhere [9]. A brief outline 
of the components that are most relevant for this analysis is given below. 



The charge-conjugated modes are implied throughout this paper. 
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In the kinematic range of the analysis, charged particles were tracked in the central track- 
ing detector (CTD) [10] and the MVD [11]. These components operated in a magnetic 
field of 1.43 T provided by a thin superconducting solenoid. The CTD consisted of 72 
cylindrical drift chamber layers, organised in nine superlayers covering the polar-angl^ 
region 15° < ^ < 164°. 

The MVD consisted of a barrel (BMVD) and a forward (FMVD) section with three 
cylindrical layers and four planar layers of single-sided silicon strip sensors in the BMVD 
and FMVD respectively. The BMVD provided polar-angle coverage for tracks with three 
measurements from 30° to 150°. The FMVD extended the polar-angle coverage in the 
forward region to 7°. After alignment, the single- hit resolution of the BMVD was 25 /xm 
and the impact-parameter resolution of the CTD-BMVD system for high-momentum 
tracks was 100 fim. 

The high-resolution uranium-scintillator calorimeter (CAL) [12] consisted of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part was 
subdivided transversely into towers and longitudinally into one electromagnetic section 
(EMC) and either one (in RCAL) or two (BCAL and FCAL) hadronic sections (HAC). 
The smallest subdivision of the calorimeter was called a cell. The CAL energy resolutions, 
as measured under test-beam conditions, were a{E)/E = 0.18/v^ for electrons and 
a{E)/E = 0.35/ VE for hadrons, with E in GeV. 

The position of the scattered electron was determined by combining information from the 
CAL and, where available, the small-angle rear tracking detector (SRTD) [13] and the 
hadron-electron separator (HES) [14]. 

The luminosity was measured using the Bethe-Heitler reaction ep — e'jp with the lu- 
minosity detector which consisted of independent lead-scintillator calorimeter [15] and 
magnetic spectrometer [16] systems. The fractional systematic uncertainty on the mea- 
sured luminosity was 2.6%. 



3 Event selection and reconstruction 

A three-level trigger system was used to select events online [9, 17, 18]. At the third level, 
events with a reconstructed scattered-electron or D-meson candidate were kept for further 
analysis. 



^ The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. 
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The kinematic variables Q^, x and the fraction of the electron energy transferred to the 
proton in its rest frame, y, were reconstructed using the double angle (DA) method [19] 
which relies on the angles of the scattered electron and the hadronic energy flow. 

The events were selected offline with the following cuts: 

• Ee' > lOGeV, where E^' is the energy of the scattered electron; 

• He < 0.95, where Ue is determined from the energy and angle of the scattered electron. 
This condition removes events where fake electrons are found in the FCAL; 

• 2/JB > 0.02, where JB signifies the Jacquet-Blondel [20] method of kinematic recon- 
struction. This condition rejects events where the hadronic system cannot be measured 



• 40 < 5 < 65 GeV, where 5 = ^ — cos(6'j)) and Ei is the energy of the i^^ energy- 
fiow object (EFO) [21] reconstructed from tracks detected in the CTD and MVD and 
energy clusters measured in the CAL. The sum i runs over all EFOs; 

• l^vtxl < 50 cm, where Z^tx is the primary vertex position determined from tracks; 

• the impact point (X, y), of the scattered electron on the surface of the RCAL must 
lie outside the region (±15 cm, ±15 cm) centred on (0,0). 

Electron candidates in the transition regions between FCAL and BCAL as well as between 
BCAL and RCAL were rejected because of the poor energy reconstruction in these areas. 
The angle of the scattered electron was determined using either its impact position on the 
CAL inner face or a reconstructed track. When available, SRTD and HES were also used. 
The energy of the scattered electron was corrected for non-uniformity due to geometric 
effects caused by cell and module boundaries. 

The selected kinematic region was 5 < < 1000 GeV^ and 0.02 < y < 0.7. The pro- 
duction of and mesons was measured in the range of transverse momentum 
1.5 < < 15 GeV and pseudorapidity \ri^\ < 1.6. 

The decay-length significance is a powerful variable for the rejection of combinatorial 
background and is defined a.s Si = l/ai, where I is the decay length in the transverse plane 
and ai is the uncertainty associated with this distance. The decay length is the distance 
in the transverse plane between the point of creation and decay vertex of the meson and 
is given by 



where pf. is the transverse momentum vector and Sxy is the two dimensional position 



precisely; 




(1) 



vector of the reconstructed decay vertex projected onto the XY plane. The vector Bxy 
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points to the fitted geometrical centre of tlie beam-spot wliicli is taken as tlie origin of tlie 
D meson. Tlie centre of the elliptical beam-spot was determined every 2000 well measured 
events [22] by fitting a Gaussian curve to the X, Y and Z distributions of the primary 
vertex. The mean of these fitted curves was then taken to be the beam-spot position. The 
widths of the beam-spot were 80 /xm and 20 /im in the X and Y directions, respectively. 
The decay- length error, o"/, was determined by folding the width of the beam-spot with 
the covariance matrix of the decay vertex after both were projected onto the D meson 
momentum vector. 

3.1 Z)-meson reconstruction 

The (and D~) mesons were reconstructed in the decay channel /i ~7r+7r+(+c.c.). 

In each event, all track pairs with equal charges were combined with a third track with 
opposite charge to form a candidate. The pion mass was assigned to the tracks with 
equal charges and the kaon mass was assigned to the remaining track. These were then 
associated and refitted to a common decay vertex [23] and the invariant mass, MlK-Kir), 
was calculated. The tracks were required to have transverse momentum p'!^ > 0.25 GeV 
and pff > 0.5 GeV for the pion and kaon tracks, respectively. To ensure that all tracks 
used were well reconstructed they were required to have passed through 3 superlayers 
of the CTD and have at least 2 BMVD measurements in the XY plane and 2 in the Z 
direction. 

Figure [T] shows the M{K'ktt) distribution for candidates. The combinatorial back- 
ground was reduced by the requirements that the of the decay vertex be less than 9 for 
3 degrees of freedom and that the decay-length significance, S"/, be greater than 3 (see Fig. 
H]). In order to extract the number of reconstructed mesons the M{K7nr) distribution 
was fitted with the sum of a modified Gaussian function [24] and a linear background 
function. The modified Gaussian function used was 

Gauss'"°^ oc exp [-0.5- a;i+i/(i+o.5.x)j ^ ^2) 

where x = | [M(i^ vrvr) — Mo]/(t|. This functional form described both the data and MC 
well. The signal position, Mq, and the width, a, as well as the numbers of mesons 
in each signal were free parameters of the fit. The number of reconstructed mesons 
yielded by the fit was N{D+) = 3995 ± 156. 

A sample of candidates with > 0.5 GeV, pff > 0.7 GeV and > 3 GeV was 
used to obtain the lifetime of the meson. The higher p^ cuts were used to obtain a 
signal with no requirements made on the significance of the decay length. The number of 
reconstructed mesons yielded by the fit to the data was N{D~^) = 4383 ± 353. 
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The D° (and D^) mesons were reconstructed in the decay channel —>■ i^~7r+(+c.c.), 
with candidates found in a similar manner to the D^, except that only oppositely charged 
pairs of tracks were combined together to form the meson candidate. The tracks were 
required to have transverse momentum pff > 0.7 GeV and > 0.3 GeV for the kaon and 
pion tracks. The ^"^^ Si cuts were 8 and 1, respectively, with 1 degree of freedom in 
the vertex fit (see Fig. Hj). After selection, the candidates were separated into tagged 
and antitagged samples with the antitagged sample used for cross-section measurements. 

The tagged group consisted of candidates which are consistent with a D*"^ D^T^f 
decay when combined with a third track that could be a "soft" pion, (vTs). The soft 
pion was required to have px > 0.12 GeV and charge opposite to that of the kaon. The 
tagged sample was used for the correction of the MC and reflection subtraction in the 
antitagged sample. For the antitagged sample, containing mesons not coming from a 
D*"^, incorrect assignment of the pion and kaon masses produced a wider reflected signal. 
The distribution of this reflection was estimated using the tagged candidates and, 
after normalising it to the ratio of the number of mesons in the two samples, it was 
subtracted from the antitagged candidates. Figure [2] shows the M{Ktt) distributions 
for tagged and antitagged candidates. The distributions were fitted simultaneously 
assuming that both have the same peak position and width and, like the , were param- 
eterised as a modified Gaussian function. The number of antitagged (tagged) mesons 
yielded by the fit was A^^^^titagpo) ^ 5534 ^ 345 (Artagpo^ ^ ^ggg ^ ^j^y 

A sample of candidates with p^^ > 0.8 GeV and > 3 GeV was used to obtain the 
lifetime of the meson. The higher pt cuts were used to obtain a signal with no require- 
ments made on the significance of the decay length. The number of antitagged (tagged) 
L)0 mesons yielded by the fit was Ar^ntitagpo) ^ 55^2 ± 283 (A^*^s(L)0) = 1495 ± 56). 

3.2 Z)-meson lifetimes 

Lifetimes for the and mesons were calculated using decay lengths in the transverse 
plane and reconstructed D-meson signals in the kinematic region 5 < < 1000 GeV^, 
0.02 < y < 0.7, 3 < pfj^ < 15 GeV and \r]^\ < 1.6. Unfolding is not necessary as the 
detector acceptance is uniform with respect to the displacement of the secondary vertex 
and the normalisation of the lifetime distribution is irrelevant. The number of D mesons 
in a given bin of proper decay length, ct, was extracted and the distributions fitted with 
the function 

00 

j e-'^'du, (3) 

^min 



fict) 



2A 



exp 





fct 








--) 
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where Uram = {—ct/a + a/X), X is the hfetime and a is the spatial resolution. This 
function represents an exponential decay convoluted with a Gaussian resolution. For the 
purposes of the lifetime extraction, a was set to the value extracted from the tagged 
sample, 120 /im, which depended only weakly on pt- 

The fitted ct distributions for and mesons are shown in Fig. [3] and the extracted 
values for the lifetime are: 



The systematic uncertainties are significantly smaller than the statistical uncertainty 
as the measurement has only a small dependence on the details of the MC simula- 
tion. The values are consistent with the world average values of 311.8 ± 2.1 //m and 
122.9 ± 0.5 /im [25] for the D+ and D^, respectively. 

4 Monte Carlo models 

The acceptances were calculated using the Rapgap 3.00 [26] Monte Carlo (MC) model, 
which was interfaced with Heracles 4.6.1 [27] in order to incorporate first-order elec- 
troweak corrections. The generated events were then passed through a full simulation of 
the detector using Geant 3.21 [28] before being processed and selected with the same 
software as used for the data. 

The MC was used to simulate events containing charm produced by the BGF process. The 
Rapgap generator used leading-order matrix elements with leading- logarithmic parton- 
shower radiation. The CTEQ5L [29] PDF for the proton was used, and the charm quark 
mass was set to 1.5 GeV. Charm fragmentation was implemented using the Lund string 
model [30]. D mesons originating from B decays were accounted for by inclusion of a 
Rapgap 6-quark sample where the 6-quark mass was set to 4.75 GeV. 

A weighting procedure utilising the tagged sample was applied in order to correct for 
imperfections in the MC description of the decay-length uncertainty [22] . 



The NLO QCD predictions for the cc cross sections were obtained using the HVQDIS 
program [31] based on the fixed- flavour- number scheme (FFNS). In this scheme, only 
light partons {u, d, s, and g) are included in the initial-state proton as partons whose 



ct{D+) 
ct{D^) 



326±21(stat.)/im 
132 ± 7(stat.) /im 



5 NLO QCD calculations 
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X, Q'^ distributions obey the DGLAP equations [32] and the cc pair is produced via the 
BGF mechanism with NLO corrections [33]. The presence of different large scales, Q, 
Pt and the mass of the c quark, m^, can spoil the convergence of the perturbative series 
because the neglected terms of orders higher than (where as is the strong coupling 
constant) contain \og{Q'^/m'^) factors which can become large. 

The predictions for D-meson production at NLO were obtained using HVQDIS with the 
following inputs. The ZEUS-S NLO QCD global fit [34] to structure function data was 
used as the parameterisation of the proton PDFs. This fit was repeated [35] in the FFNS, 

(3) 

in which the PDF has three active quark flavours in the proton. In this fit Aq^j-, was 
set to 0.363 GeV and the mass of the charm quark was set to 1.5 GeV; the same mass 
was therefore used in the HVQDIS calculation. The renormalisation and factorisation 
scale, /i = = /iF, was set to ^/Q"^ + Ami. The charm fragmentation to the particular 
D meson was described by the Peterson function [36] with the Peterson parameter, e, 
set to 0.035 [37]. The values used for the hadronisation fractions, /(c D), were those 
previously measured in DIS at ZEUS, 0.216i'^o;q29 and 0.450l;o;Qgo for the and antitagged 
D^, respectively [8]. 

To estimate the contribution of beauty production, the HVQDIS calculation and hadro- 
nisation from the MC were combined, using da{b -D)nlo+mc = da{bb)^i,o- Chad where 
Chad = da{b D)Mc/d(j{bb)Mc- The ZEUS NLO QCD fit was used as the proton PDF, so 
that the mass used in this fit was also used in the HVQDIS program. The hadronisation 
fraction, f{b D), was set to 0.231 and 0.596 for the and D^, respectively [38]. 

The HVQDIS predictions for D-meson production are affected by theoretical uncertainties 
listed below. The average uncertainty on the total cross sections is given in parentheses: 

• the ZEUS PDF uncertainties propagated from the experimental uncertainties of the 
fitted data (±5%). The change in the cross section was independent of the kinematic 
region; 

• the mass of the charm quark (±8%). The charm quark mass was changed consistently 
in the PDF fit and in HVQDIS by ^0.15 GeV; 

• the renormalisation and factorisation scale, The scales 2^/(5^ + 4m^ and 
a/(5^/4 + ml were used; 

• the e parameter of the Peterson fragmentation function (ly^) was varied by l^ois [39]. 
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6 Data correction and systematic uncertainties 



For a given observable Y, the production cross section was determined using: 



da 
dY 



N{D) 



(4) 



A-C-B-AY 



where N{D) is the number of reconstructed D mesons in a bin of size AY, A is the 
reconstruction acceptance as found from the MC sample which includes migrations, effi- 
ciencies and QED radiative effects for that bin, £ is the integrated luminosity and B is 
the branching ratio for the decay channel used in the reconstruction. 

Small admixtures to the reconstructed signals from other decay modes were taken into 
account in the MC sample used for the acceptance-correction procedure. To correct from 
the number of reconstructed mesons to the production cross sections, small migrations 
were taken into account between the tagged and antitagged samples. It was checked that 
the Rapgap MC sample gives a reasonable description of the data for selected DIS and 
D meson variables. Figures HI [5] and [6] show important variables for the secondary vertex 
reconstruction, distributions for the DIS variables and the kinematics of the D meson, 
respectively. For all variables, the number of reconstructed D mesons is extracted by 
fitting the number of D mesons in each bin of the distribution. The MC provides a good 
enough description of the data for acceptance calculations in all variables. 

Reconstruction acceptances vary depending on the particle and kinematic region of the 
measurement. For example, the overall and D° acceptances calculated with Rapgap 
after applying the selection criteria for the kinematic region are ~ 7% and ~ 17%, re- 
spectively. The lower average acceptance in relation to previous ZEUS measurements is 
accounted for by reduced efficiency for reconstructed D mesons due to the extension of 
the kinematic range to lower p!^ and the use of lifetime tagging. This is offset by a gain 
of a factor of 20 and 3 in the signal to background ratios of the -D+ and samples. 

The systematic uncertainties of the measured cross sections were determined by changing 
the analysis procedure and repeating all calculations. The following possible sources of 
systematic uncertainties were considered [22,40] with the average effect on the measured 
and total cross sections shown in parentheses: 

• {6i} the cut on yjB was changed by ^.o^ot 

• {62} the cut on the scattered electron energy, E^', was changed by ±1 GeV (1^%); 

• {6^} the 6-quark cross section was varied by a factor of two in the reference MC sample 



• {64} the uncertainty of the tracking performance was obtained by varying all momenta 
by ±0.3% which corresponds to the uncertainty in the magnetic field; and by changing 



+1.3% V 
-1.7%^' 
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the track momentum and angular resolutions by of their values. The asymmetric 
resolution variations were used since the MC signals typically had somewhat narrower 
widths than those observed in the data (±1%); 

• {^s} the uncertainty of the MVD hit efficiency was obtained [22] by evaluating the 
relative difference in single-track efficiency between data and MC when 2 XY and 2 
Z measurements were required in the BMVD (±1.1%); 

• {5q} the cut on Si was varied by ±1.0 in the analysis and ±0.4 in the analy- 
sis [41] 

• {5-j} the cut on the of secondary vertex was changed by ±2 in the analysis 
and ±1.5 in the analysis (^^%); 

• {^s} the MC distribution was reweighted in order to account for the difference (see 
Fig. E]) between data and MC (< 1%); 

• {5q} the MC 7]^ distribution was reweighted in order to account for the difference (see 
Fig. ED between data and MC (< 1%). 

An additional source of systematic uncertainty in the analysis was investigated: 

• {^10 } the background function was parameterised by an exponential function (±4%). 

Several other sources of systematic uncertainty were considered and found to have an 
effect of < 1% on the total cross sections. These sources were related to the DIS selection 
criteria and the method for extracting the number of tagged mesons. 

The systematic uncertainty is dominated by ^e, which is related to the description of 
the MVD resolution. This uncertainty was evaluated from the differences between the 
data and MC description of ai (see Fig. H]). This difference was then propagated to a cut 
variation of the Si cut and the analysis procedure repeated. 

Contributions from the different systematic uncertainties were calculated and added in 
quadrature separately for positive and negative variations. Uncertainties due to those 
on the luminosity measurement and branching ratios were only included in the measured 
and total cross sections. For differential cross sections these uncertainties are not 
included in the tables and figures. 

7 Cross sections 

Charm meson cross sections for the process ep — > eDX were calculated using the recon- 
structed and signals (see Section [3]) in the kinematic region 5 < < 1000 GeV^, 
0.02 <y < 0.7, 1.5 < p^ < Ih GeV and |r/^| < 1.6. 

The following cross sections were measured: 
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• The production cross section for and D mesons: 

cr(D+) = 4.67 ±0.26 (stat.) (syst.) ± 0.17(br.) ± 0.12(lumi.) nb 

• The production cross section for and mesons not originating from D*^ decays: 

^antitag^^O) = 7.49 ± 0.46 (stat.) toil (syst.) ± 0.14(br.) ± 0.20(lumi.) nb 

The corresponding predictions from HVQDIS are: 

= 4.42 1°:^^ (syst.) t^:^^ (had.) nb 

^antitagpO) ^ g 35 +1-79 (gyg^ ) +0.52 ^^^^ -^ 

where "had." and "br." represent the uncertainty on the HVQDIS prediction due to the 
uncertainties of the hadronisation fraction /(c D) and decay-chain branching ratios, 
respectively. The predictions used the defauh parameter settings as discussed in Section 
m The quadratic sum of the other uncertainties of these predictions is shown with the 
"syst." label. A small contribution (~ 2%) to the total cross sections arises from D 
mesons produced in bb events. All predictions include a bb contribution calculated in each 
bin with HVQDIS. The HVQDIS predictions are in agreement with the data. 

The differential and cross sections as functions of Q'^,x,pi^ and rj^ are shown in 
Figs. [7] and M and given in Tables [T] and [H The cross sections in and x both fall by 
about three orders of magnitude in the measured region. The cross section in falls 
by about two orders of magnitude and there is no significant dependence on t]^. The 
HVQDIS predictions describe the shape of all measured differential cross sections well. 
The slight difference in normalisation in Fig. [8] reflects the difference of the corresponding 
total cross section. 



8 Extraction of F^^ 

The open-charm contribution, F^*^, to the proton structure function, F2, can be defined 
in terms of the inclusive double-differential cc cross section in x and by 

"'^l^^f = 5^ { [1 + (1 - Ff(x, Q') - y^F-ix, Q')} . (5) 

In this paper, the cc cross section is obtained by measuring the and production 
cross sections and employing the hadronisation fraction /(c D) to derive the total 
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charm cross section. A limited kinematic region is accessible for the measurement of D 
mesons, therefore a prescription for extrapolating to the full kinematic phase space is 
needed. The measured value of Fl'^ in a bin i is calculated with 



■^2,racas{^i^Qi) ~ _ / , n -^2,thco ) ) (6) 



^2^ _ ^i,meas(ep ~^ DX) 

o-i,theo(ep DX) 

where cij^meas is the cross section in the bin i in the measured region of pf? and rj^ and 
c"i,theo is the corresponding cross section evaluated with HVQDIS. The value of i^l^theo 
was calculated in FFNS from the NLO coefficient functions [34] using the same values 
of parameters as in the calculation of Cj^thco- The cross sections cTj j^eas(ep DX) were 
measured in bins of and y (Table [3]) and Fg'^ is quoted at representative and x 
values near the centre-of-gravity for each bin (Table H]). 

Beauty contributions were subtracted from the data using the predictions obtained from 
HVQDIS. The contribution to the total cross section from Ff^ calculated using the ZEUS 
NLO fit was, on average, 1.3% and at most 4.7% [5] and was taken into account, in CTj^theo, 
in the extraction of F2^. The size of the contribution from was similar to that in other 
PDFs. 

The factor to extrapolate from the measurement range to the full phase space was esti- 
mated using HVQDIS and was found to vary from f» 1.5 at high to ~ 3.2 at low Q"^. 
A complete list of the extrapolation factors is given in Table ID 

The following uncertainties associated with the method of extrapolation were evaluated 
with the average effect given in parentheses: 

• changing the charm mass by =F0.15GeV consistently in the HVQDIS calculation and 
in the calculation of F2'\yico (±2%). The largest effect was seen at low x and low 

• using the upper and lower predictions given by the uncertainty in the ZEUS NLO PDF 
fit, propagated from the experimental uncertainties of the fitted data, to perform the 
extraction of Ff (< 1%); 

• changing the contribution of beauty events subtracted from the data by a factor 2 
i1ll%). The largest effect was seen at low x and high (^7%); 

• a Lund string model in Rapgap was used as in previous analyses [1,5,8] rather than 
the Peterson function in HVQDIS (±7%). The largest effect was seen at high x and 
low g2 (±14%). 

The ^2^^ values measured from D~^ and D'^ data are combined using a procedure that 
accounts for the systematic and point-to-point correlations between the analyses [42]. 
The combined values of F2'^ obtained from D~^ and D^ production are given in Table [5] 
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and shown in Fig. O Also shown is the ZEUS NLO QCD fit which describes the data 
well for all and x. The uncertainty of the theoretical prediction shown is that from 
uncertainty of the charm mass. Due to the improved statistical precision resulting from 
lifetime tags with the MVD, more measurements of were extracted than in the previous 
publication [8]. Also, extrapolation factors were significantly reduced, from e.g. a value 
of about 5-6 to about 2 at = 20.4 GeV^, due to the extension of the kinematic range 
to lower pJ^. At high these results are competitive with Z)*^ based measurements [5]. 

9 Conclusions 

The production of the charm mesons and has been measured with the ZEUS 
detector in the kinematic range 5 < < 1000 GeV^, 0.02 <y < 0.7, 1.5 < < 15GeV 
and |?7^| < 1.6. Combinatorial background to the D meson signals was reduced by using 
the ZEUS microvertex detector to reconstruct displaced secondary vertices. 

The measured D meson cross sections were compared to the predictions of NLO QCD 
with the proton PDFs extracted from inclusive DIS data. A good description was found. 

The visible cross sections in bins of y and were used to extract the open-charm contri- 
bution, Fl'^, to the proton structure function, F2. The extraction used factors calculated 
within the framework of NLO QCD. 

The use of the microvertex detector has increased the precision and allowed an extension 
in the kinematic range to lower values of compared to previous results. Along with 
previous measurements of F^^^, the results presented here provide a direct constraint on 
the gluon density of the proton. 
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Table 1: Measured cross sections as a function of Q^, x, p^ and rj^ for 
5 < g2 < 1000 GeV^, 0.02 <y < 0.7, 1.5 < < 15 GeV and |r/^*| < 1.6. The 
statistical and systematic uncertainties are shown separately. The cross sections 
have further uncertainties o/3.5% from the Ji ^7r+7r+(±c.c.) branching ratio, 

and 2.6% from the uncertainty in the luminosity measurement. 
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Table 2: Measured cross sections for D^/ not coming from a D*^ as a function 

ofQ^, X, and r]^'/^" for 5 < < 1000 GeV^, 0.02 < y < 0.7, 1.5 < 

p!^ < 15 GeV and \ri^°^^°\ < 1.6. The statistical and systematic uncertainties 
are shown separately. The cross sections have further uncertainties of 1.9% from 
the — > K~'K^{+c.c.) branching ratio, and 2.6% from the uncertainty in the 
luminosity measurement. 
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0.20 ±0.05 +0.01 -0.03 
0.35 ±0.06 +0.05 -0.08 
0.24 ±0.05 +0.03 -0.06 




g2 bin 
(GeV^) 


y bin 


(nb) 


5, 9 


0.02, 0.12 
0.12, 0.30 
0.30, 0.70 


0.80 ±0.24 +0.23 -0.16 
0.92 ±0.20 +0.13 -0.12 
0.48 ±0.17 +0.11 -0.14 


9, 44 


0.02, 0.12 
0.12, 0.30 
0.30, 0.70 


1.62 ±0.18 +0.10 -0.13 
1.42 ±0.15 +0.05 -0.06 
1.22 ±0.24 +0.24 -0.18 


44, 1000 


0.02, 0.12 
0.12, 0.30 
0.30, 0.70 


0.19 ±0.09 +0.06 -0.04 
0.54 ±0.09 +0.06 -0.04 
0.54 ±0.15 +0.15 -0.18 



Table 3: Measured cross sections for and D^/D^ not coming from a D*^ 
m each of the and y hms for 5 < < iqoo GeV^ , 0.02 < y < 0.7, 1.5 < 
^p^ < 15 GeV and \ri^\ < 1.6. The statistical and systematic uncertainties are 
shown separately. The and D^/D^ cross sections have further uncertainties of 
3.5% and 1.9% from the — > K~t[^t[^[+c.c.) and — K~tt^{+c.c.) branching 
ratios. The additional uncertainty from the luminosity measurements is 2.6%. 
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(GeV^) 


X 


^stat ^syst 


^cxtrap factor 


7.0 


0.00022 
0.00046 
0.00202 


0.295 ±0.092 +0.091 -0.074 
0.176 ±0.031 +0.023 -0.050 
0.091 ±0.023 +0.030 -0.025 


+0.026 -0.022 3.2 
+0.010 -0.008 2.3 
+0.013 -0.014 3.1 


20.4 


0.00065 
0.00134 
0.00588 


0.319 ±0.054 +0.037 -0.086 
0.241 ±0.024 +0.016 -0.014 
0.131 ±0.015 +0.010 -0.018 


+0.022 -0.020 2.5 
+0.013 -0.013 1.8 
+0.009 -0.009 2.4 


112.0 


0.00356 
0.00738 
0.03230 


0.260 ±0.058 +0.029 -0.066 
0.280 ±0.049 +0.038 -0.064 
0.089 ±0.024 +0.004 -0.015 


+0.020 -0.025 1.7 
+0.032 -0.033 1.5 
+0.002 -0.002 2.4 



(GeV2) 


X 


Agtat Asyst 


^extrap factor 


7.0 


0.00022 
0.00046 
0.00202 


0.116 ±0.042 +0.028 -0.035 
0.131 ±0.029 +0.019 -0.017 
0.068 ±0.020 +0.019 -0.014 


+0.010 -0.009 3.2 
+0.007 -0.006 2.3 
+0.010 -0.010 3.1 


20.4 


0.00065 
0.00134 
0.00588 


0.252 ±0.051 +0.049 -0.037 
0.169 ±0.019 +0.006 -0.007 
0.109 ±0.012 +0.006 -0.009 


+0.017 -0.016 2.5 
+0.009 -0.009 1.8 
+0.007 -0.008 2.4 


112.0 


0.00356 
0.00738 
0.03230 


0.280 ±0.086 +0.077 -0.096 
0.203 ±0.037 +0.024 -0.016 
0.040 ±0.019 +0.012 -0.008 


+0.022 -0.027 1.7 
+0.023 -0.024 1.5 
+0.001 -0.001 2.4 



Table 4: The extracted values of from the production cross sections of D and 
D^/D^ not coming from D*^ at each and x value. The statistical, systematic 
and extrapolation uncertainties are shown separately. The values of the extrapola- 
tion factor used to correct to the full p^ and rj^ phase space are also shown. The 
values extracted from and D^/D^ have further uncertainties as detailed in the 
caption to Tahle\^ 
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(GeV2) 


X 


F2 Astat Asyst 


Aextrap 


7.0 


0.00022 
0.00046 
0.00202 


0.260 ±0.062 ±0.091 
0.157 ±0.022 ±0.031 
0.088 ±0.017 ±0.028 


±0.007 -0.067 
±0.016 -0.035 
±0.009 -0.016 


20.4 


0.00065 
0.00134 
0.00588 


0.291 ±0.038 ±0.064 
0.213 ±0.016 ±0.014 
0.126 ±0.010 ±0.014 


±0.020 -0.094 
±0.018 -0.040 
±0.010 -0.042 


112.0 


0.00356 
0.00738 
0.03230 


0.257 ±0.046 ±0.057 
0.238 ±0.030 ±0.039 
0.086 ±0.020 ±0.018 


±0.020 -0.084 
±0.015 -0.041 
±0.001 -0.026 



Table 5: The combined Fg'^ values from the production cross sections of D 
and D^/D^ not coming from D*^ at each and x value. The statistical, 
systematic and extrapolation uncertainties are shown separately. The measure- 
ments have a further uncertainty of 3.3% from the — »• K~7i~^7t~^{+c.c.) and 
— >• K~Ti^{+c.c.) branching ratios. The additional uncertainty from the lumi- 
nosity measurement is 2.6%. 
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Figure 1: The M^Knn) distribution for the candidates (dots). The solid 
curve represents a fit to the sum of a modified Gaussian function and a linear 
background function. 
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Figure 2: The M{Ktt) distributions (dots) for (a) D^/D^ candidates not con- 
sistent with a D*^ decay, obtained after the reflection subtraction (see text) and 
(b) D^/D^ candidates consistent with a D*^ decay. The solid curves represent the 
simultaneous fit as described in the text. 
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Figure 3: The distributions of reconstructed candidates (circles) andD^/D^ 
candidates not consistent with a D*^ decay (triangles) extracted in bins of proper 
decay length, ct. Both distributions are fitted with functions described by a Gaussian 
convoluted with an exponential decay. The relative normalisation of the distribu- 
tions is chosen to aid visibility. 
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Figure 4: Reconstructed decay-vertex variables Si and for (a, b) and (c, d) 
D^/D^ mesons. Data (points) are compared to detector-levelllAPGAP predictions 
(shaded histograms). All histograms are normalised to unit area. The dashed line 
indicates regions removed by the cuts placed on these variables. 
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Figure 5: Reconstructed DIS variables for events with candidates (extracted 
with a fitted signal) for data (points) compared to detector-level Rapgap predictions 
(shaded histograms) . All histograms are normalised to unit area. Similar agreement 
was observed for the D^/D^ candidates. 
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Figure 6: Reconstructed (a, b) and (c, d) D^/D^ kinematic variables for 
data (points) compared to detector-level Rapgap predictions (shaded histograms). 
All histograms are normalised to unit area. 
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Figure 7: Differential cross sections for mesons as a function of (a) , (b) 
X, (c) and (d) r]^ compared to the NLO QCD predictions of HVQDIS. Sta- 
tistical uncertainties are shown by the inner error bars. Statistical and systematic 
uncertainties added in quadrature are shown by the outer error bars with the shaded 
region representing the uncertainty of the HVQDIS prediction. The ratios, R, of 
the cross sections to the central HVQDIS prediction are also shown in the lower 
section of each plot. 
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Figure 8: Differential cross sections for D^/D^ mesons not from D*^ decay 



as a function of (a) Q , (h) x, (c) pj. 



DO /DO 



and (d) r]^"/^" compared to the NLO 



QCD predictions of HVQDIS. Statistical uncertainties are shown by the inner error 
bars. Statistical and systematic uncertainties added in quadrature are shown by the 
outer error bars with the shaded region representing the uncertainty of the HVQDIS 
prediction. The ratios, R, of the cross sections to the central HVQDIS prediction 
are also shown in the lower section of each plot. 
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Figure 9: Combined values 0/^2^*^ extracted from and / not from D*^ 
(circles) as a function of x in three bins of Q"^. The data are shown with statis- 
tical uncertainties (inner bars) and statistical and systematic uncertainties added 
in quadrature (outer bars) and, where possible, are compared to previous ZEUS 
measurements with these mesons. The measurements have a further uncertainty of 
3.3% from the — > K'tt^tt^^+c.c.) and K~'k^{+c.c.) branching ratios. 

The additional uncertainty from the luminosity mesaurements is 2.6%. The shaded 
band shows the predicted values of F2'^ for values ofrric between 1.35 and 1.65 GeV. 
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